A Cautionary Note on Cosmological Magnetic Fields 
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This note is concerned with potentially misleading concepts in the treatment of cosmological 
magnetic fields by magnetohydrodynamical (MIfD) modelling, ft is not a criticism of MHD itself 
but rather a cautionary comment on the validity of its use in cosmology. Now that cosmological 
data are greatly improved compared with a few decades ago, and even better data are imminent, it 
makes sense to revisit original modelling assumptions and examine critically their shortcomings in 
respect of modern science. Specifically this article argues that ideal MHD is a poor approximation 
around recombination, since it inherently restricts evolutionary timescales, and is often misapplied 
in the existing literature. 

PACS numbers: 98.80.-k, 98.80.Jk, 91.25.Cw, 94.20.wc, 94.30.cs, 95.30.Qd 



I. INTRODUCTION 

The role of the magnetic field in the evolution of 
the Universe has been a challenging problem for many 
decades, retaining a frustrating degree of speculation de- 
spite remarkably inventive and ingenious mathematical 
modelling. Definitive progress to date has been ham- 
pered by the ambivalence of observational data: as yet, 
contradictory hypotheses continue to remain possibilities 
for as long as the definitive observation remains elusive. 

However, the quantity and precision of appropriate 
cosmological data are soon to be revolutionised with the 
imminent construction of dedicated instruments such as 
SKA, and the era of precision diagnosis of cosmological 
magnetic fields is not far away. It is surely appropri- 
ate then to re-examine the theoretical basis upon which 
the modelling of cosmological magnetic fields is based, in 
preparation for these new data. 

A favourite mathematical framework for modelling the 
evolution of magnetism near the decoupling era (and 
sometimes for the early universe) is magnetohydrody- 
namics (MHD) , which blends the continuum properties of 
the matter in the Universe with a restricted set of electro- 
magnetic concepts, yielding a hybrid magneto-fluid con- 
text in which to study the influence of matter on mag- 
netic fields, and vice- versa. Such MHD models are often 
imbued with perfect electrical conductivity (that is, zero 
electrical resistivity), and are almost invariably single- 
fiuid, since this combination allows significant simplifica- 
tion of mathematical modelling, without apparently sac- 
rificing too much of the underlying physics. It has to be 
mentioned that many of the cosmological MHD models 
used in this context have additional properties not en- 
countered in standard plasma modelling; these primarily 
arise from the incorporation of Universal expansion, and 
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sometimes from accommodating radiation pressure as an 
extra force term. 

In this short article we examine critically the physical 
underpinning of the magnetofiuid Universe as a suitable 
model for cosmological magnetogenesis and evolution, 
and find such models to be flawed in basic concept. We 
conclude that MHD, and its variants, are quite unsuited 
to modelling the critical phenomena; instead we believe 
that consistency demands that we use a relativistic ki- 
netic theory instead [see fll| for a gas kinetics (but not 
plasma) treatment in a Friedman-Lemaitre-Robertson- 
Walker (FLRW) space-time]. The next sections explain 
our reasoning, by highlighting the inadequacies of the 
magnetofiuid approach, ranging over issues such as the 
calculation of conductivity, the meaning of the current 
density, the interpretation of mass density and the allied 
velocity field, and finally comments on the juxtaposition 
of relativity and electromagnetism with an MHD plasma 
description. 



II. THE RELEVANCE OF MHD TO 
COSMOLOGY 

Magnetohydrodynamics is a fiuid model, in which as- 
pects of electromagnetism are incorporated into standard 
hydrodynamics to give an electrically conducting contin- 
uum that can generate a self-field, and respond to an 
applied one. Single fluid MHD cannot support charge 
separation, and perfectly conducting (or ideal) MHD has 
zero electric field in the rest frame of the plasma. 

MHD is the favoured mathematical model for the de- 
scription of the cosmic plasma behaviour but caution is 
required in the interpretation of the results. MHD is an 
excellent plasma model in many disparate contexts, but 
care has to be taken to ensure that conclusions are not 
drawn that cannot be sustained by the restricted physics 
implicit in MHD, particularly where the context requires 
elements of relativity. 



In cosmological MHD there may be at best a conflict in 
terminology, and at worst, flawed physics. In order to be 
clear about the problems, we elaborate below on general 
issues in which the physical interpretation is perhaps at 
odds with the mathematical framework. 

One fundamental aspect in cosmological descriptions 
involving plasma is as follows: in the literature it is 
assumed that all the universe is ionised, meaning the 
plasma density and baryon density are identical. This 
may be reasonable before recombination (assuming that 
dark matter can be ignored) but not after! Whilst this 
might appear at first a trivial statement, the implications 
can be significant if care is not taken to identify how the 
mass densities of constituent fluids contribute to the rel- 
evant physical forces. 

Below we present a list of issues in which the use of 
MHD may lead to inconsistencies in the interpretation of 
the results. These are not in any particular order, and 
we don't claim that it is a exhaustive list. However, each 
topic highlights a specific conceptual shortcoming asso- 
ciated with the use of MHD in cosmological descriptions 
currently in the literature. 

A. The Perfectly Conducting Assumption 

One of the great strengths of single-fluid MHD is the 
simplicity in the field evolution, particularly given the 
ideal case in which the plasma has zero resistivity. Given 
that the single-fluid model immediately rules out charge 
separation, the further imposition of a perfectly conduct- 
ing fluid allows a very simple Ohm's law to be placed at 
the heart of the physics: 

E + uxB = (1) 

where E, u and B are respectively the plasma elec- 
tric field, bulk fluid velocity and magnetic flux density 
(loosely, magnetic field). One consequence is immediate: 
the electric field can only be produced as the result of a 
frame change, and is therefore entirely prescribed by the 
velocity and magnetic fields, having no independent evo- 
lution. (There are no induced fields, for example.) This 
is because any applied electric field invokes a plasma dy- 
namical response in terms of the particle current density 
J consistent with the relation 

J = crE (2) 

where a is the conductivity tensor; not a scalar, since the 
magnetic field changes the transport properties of mobile 
species for cross-field motion compared to the isotropic 
case. Clearly then unlimited particle currents can be 
sustained in the perfect conductivity case by vanishingly 
small electric fields; since such currents are not seen, then 
the Universal electric field must be zero except for that 
arising from a frame translation of the magnetic field. 
Whilst this is perfectly acceptable within the classical 
concept of MHD plasmas, it pays to consider carefully 
the implications of this in a cosmological context. 



In the early Universe (at least pre-decoupling) , the 
plasma (charged baryons and electrons) co-exists with 
a radiation field that interacts with it via the Thompson 
scattering of photons by the free electrons (to identify but 
one process). Viewed classically (that is, non quantum- 
mechanically) in terms of the propagation of the electric 
field disturbances, the plasma is a dielectric, since inter- 
action between the two continua produces a different re- 
sponse than if the radiation simply propagated in a clas- 
sical vacuum: the plasma has a refractive index different 
from unity, the classical vacuum value. A non- vacuum re- 
fractive index implies a finite conductivity. This can be 
seen from cold plasmas, for example, in which the dielec- 
tric tensor K and conductivity tensor g are intimately 
related @, [| 

K = I+—(j (3) 

where / is the unit tensor, lo is the frequency of the elec- 
tromagnetic radiation, and cq is the vacuum permittivity. 
Note that the kinetic plasma treatment is similar. Hence 
if the plasma refractive index is not unity (the classical 
vacuum value) , then the conductivity of the plasma must 
be finite, and anisotropic: particle transport parallel to 
the magnetic field is different from that perpendicular to 
it. 

Of course, in the very early universe where current car- 
riers are effectively massless then the conductivity could 
be very high [1, [^; however, the applicability of MHD in 
such energetic conditions is very much open to question. 

The transition from perfect to finite conductivity is 
subtle, but critical: finite conductivity brings evolution- 
ary and topological constraints that are simply missing 
from perfect conductivity; moreover a resistive plasma 
has additional interdependence of the fields compared 
with an ideal one. 

Note that in the recombination era, the conductiv- 
ity must evolve as the charged particle density changes, 
making any assumption of perpetual perfect conductivity 
even less credible (see the discussion later in this section). 

Hence a more appropriate general form of Ohm's Law 
is the simplest resistive one, 

E + u X B = ryJ (4) 

where rj = a^^ is the plasma resistivity, also a tensor in 
a magnetised plasma. 

It is a cosmological convention to calculate electrical 
conductivity in terms of the effect of photon scattering 
on electron transport, using a simple scalar expression 
for the conductivity based on a single interaction time t 
%]■■ 



in which r = l/(n^(Ty), with is the photon number 
density and gt is the Thomson scattering cross-section, 
and Ue is the /ree-electron number density. This yields 
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a post-recombination conductivity that is approximately 
10"'^ times that for the pre-recombination case p]. Tak- 
ing the plasma perspective, the ratio of electrical conduc- 
tivities before and after recombination is based on the 
Spitzer electrical conductivity for a fully ionized plasma 
at a temperature of 10 eV [1], and the scattering of 0.1 eV 
electrons by neutral hydrogen [l^l ■ This ratio is 



O-prc/CTpost ~ O.Oie/Xe 



(6) 



where Xe is the ionization fraction. This yields a similar 
drop in conductivity as in the Thomson scattering case if 
Xe ~ 10^^[ll|. Since the Spitzer conductivity is similar in 
value to that calculated by considering only the Thomson 
scattering of photons by electrons then this plasma cal- 
culation shows that the pure plasma treatment is equally 
valid, though the nature of the BGK approximation in 
calculating the resistivity Eq. ^ might not be wholly 
valid in either context, since it is unlikely that a single 
relaxation time can account for all scattering processes. 
Either way, the resistivity is finite on either side of recom- 
bination, and significantly lower afterwards, underlining 
the dangers in a perfectly conducting description. 



B. The Total Current 

In the present article we use the 1-1-3 covariant for- 
malism introduced by [l^, [l^l in which a general class of 
homogeneous space-times are considered. The usual V 
and t operators are generalised appropriately. 

One of the crucial approximations in MHD is that the 
current density is given exclusively in terms of the mag- 
netic field curvature: 



V X B = /ioJ 



(7) 



Note that J is not related directly to the bulk flow veloc- 
ity u, since the latter is a mass- weighted averaged quan- 
tity that is independent of charge. Hence in MHD, the 
current density is a dependent variable, and can be elim- 
inated everywhere in favour of the magnetic field. 
However, full electromagnetism demands that 



V X B 



c 



1 aE 



(8) 



where c is the speed of light, and t is the time. The ad- 
ditional term on the right is the displacement current, 
omitted from MHD because the rapid time evolution of 
electromagnetic effects is not incorporated in such a fluid 
description: there are other comparable time-dependent 
terms in the accornpanying plasma equations that have 
also been omittedj^] and so, for mathematical consis- 
tency, the displacement current must be dropped. 

The true significance of adopting Eq. ([7]) instead of 
Eq. ([8]) lies in the timcscalc for the evolution of the mag- 
netic field. With Eq. ([7]) we have 
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-GB 



(9) 



V X (u X B - 77J) - -OB 



(10) 



where we have used Eq. ^ for the most general case. Q 
denotes 3d/a, where a is the expansion factor. Eliminat- 
ing the current density using Eq. ([7|) yields 
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-GB 



(11) 



assuming the simplest case of constant resistivity. This 
parabolic equation shows that the magnetic field evolves 
on the resistive diffusion time-scale (notwithstanding the 
feedback term from the dynamo contribution, the first 
term on the right-hand side) ; cosmological expansion will 
allow an additional cosmological time-scale. 

However, if J is eliminated using Eq. ([8]) instead, the 
governing equation for the magnetic field is higher-order, 
and hyperbolic: 



Tj a^B dB 



V 



— = V X (u X B) + -^V"B - -OB (12) 

Ot .. o 



This introduces new time-scales for field evolution, in- 
cluding wave propagation, damping and growth, to be 
balanced against the dynamo term and any cosmological 
input. 

Granted, magnetic and electric fields that so arise may 
be rapidly varying: this means that the overall math- 
ematical framework must be able to accommodate fast 
time-scale evolution, a point to which we will return in a 
different context later in this article. 

Note that in some very early work (e.g. [l3|) the mag- 
netic structure was modelled without any intrinsic time 
evolution, and only varied on cosmological time scales; 
other more recent work omits the spatial structure of the 
magnetic field (e.g. 0). Neither of these approaches is 
correct: the true picture must have the correct balance 
between spatial and temporal derivatives in order to re- 
flect the correct physics. 



C. Density considerations 

One serious concern with unified fluid descriptions of 
the universe around recombination is the concept of the 
mass density. Care must be taken to distinguish be- 
tween the density of charged and non-charged particles; 
the latter are also divided between baryonic matter and 
dark matter. Only the density of plasma (that is, the 
electrically conducting and magnetised fluid) can appear 
in those plasma equations that are concerned with ex- 
clusively plasma effects. For example, neutral matter 
density cannot contribute to current densities. There 
are several examples in the literature [l^, [l^ [13] where 
only baryonic matter density enters in the model equa- 
tions with no distinction between neutrals and current- 
carrying species. Where baryonic matter is predomi- 
nantly charged, so that the baryonic and plasma mat- 
ter densities are approximately the same, this is fine. 



3 



but where neutral matter (e.g. hydrogen atoms and 
molecules) co-exists with plasma, this clearly cannot be 
correct, and care must be taken to distinguish between 
species. This latter case leads for example to an Alfven 
speed, CA = -S/^MoPOj that incorrectly uses the total 
mass density, instead of just the plasma mass density; 
the magnetic field cannot be influenced directly by the 
motion of neutral matter and vice-versa. This becomes 
particularly important when the Alfven speed is used to 
estimate magnetic field strengths in the post recombina- 
tion era (for example, pUS]). 

Momentum exchange between the fluids is obviously 
physically correct, and is a primary method of causing 
indirect interaction between the magnetic field and neu- 
tral particles, but is predicated on making the distinction 
between the different fluid types in the model [H, [l^ . 

In particular, the pressures of neutral matter and 
plasma evolve under different physical conditions. Hence, 
incorporating the magnetic field into the Jeans length re- 
quires a careful definition of the concept of magnetic and 
kinetic pressure. The latter is influenced by both neutral 
matter and plasma; the former is only directly supplied 
by the plasma. To get a true holistic picture, the neu- 
tral gas and plasma components must be identified as 
separate species from the outset, with the latter obeying 
additional force terms, but with each potentially inter- 
acting through gravitation, or explicit coupling terms. 



D. The Global Velocity Field 

This problem is related to the density issue above. 
Since the plasma and the neutral components respond 
to different dynamical equations (notwithstanding cou- 
pling) then clearly the velocity fields are different. In 
principle the neutrals need not have the same velocity as 
the plasma since they do not react directly to magnetic 
forces. Moreover, the bulk motion of neutrals makes no 
direct contribution to the evolution of magnetic field. Of 
course, should the neutral gas and the plasma be cou- 
pled in some way, then the evolution of one will affect 
the other, but this is somewhat different from assuming 
that the plasma and neutral components are somehow 
locked together. 

In the light of comments about the densities and veloc- 
ities it is appropriate to mention how plasma and non- 
plasma can be coupled to give an overall collective re- 
sponse. Plasmas and neutrals can be coupled through 
momentum transfer (Tsl . IT^ , in which each fluid exerts a 
drag on the other by virtue of relative motion. lonisation 
and recombination can also be considered as coupling 
mechanisms, in that species are converted from neutral to 
plasma and vice-versa. However single fluid MHD is not 
a good model for such processes, since charge separation 
is a basic pre-requisite but is impossible in single-fluid 
MHD plasma. Also, dark matter, neutrals and plasma 
move under the common self- gravitating potential pre- 
sented by their respective mass densities, leading to dy- 



namical equations for each species that are separate but 
coupled, emphasising that care is required in identifying 
what is meant by velocity at a given point in space since 
we have to distinguish between different species. 



E. Relativity, MHD and Photons 

In many cosmological contexts relativity is a key ele- 
ment of the physical model, including the plasma. How- 
ever, there is a problem here if the plasma is described 
by a standard fluid MHD model. Given that MHD is ex- 
clusively concerned with low frequency, long wavelength 
phenomena that are not electromagnetic in nature, it is 
far from obvious how an MHD plasma description can be 
incorporated into a Lorentz invariant model of the whole 
ensemble: MHD cannot be Lorentz invariant, since the 
displacement current has been omitted, and so the cur- 
rent density in MHD is not a Lorentz invariant. 

Furthermore, photons and free electrons are formally 
irreconcilable with an MHD prescription since MHD is 
not fully electromagnetic and is single fluid: electron and 
ions are combined in an averaged description. Two-fluid 
MHD is a possible approach but requires, amongst other 
things, a more sophisticated Ohm's law than appears in 
the literature. (Recall that single fluid MHD cannot sus- 
tain charge separation.) Hence, photon scattering using 
Thomson cross-sections cannot be rigorously quantified 
in a single MHD context since the electron number den- 
sity cannot formally be deduced independently, and the 
physics of the scattering processes involving photons and 
free electrons is not consistent with the exclusively long 
time-scale processes that are valid in MHD. Any formal 
attempt to combine both approaches in a unified treat- 
ment of the damping of MHD waves (for example in [I7|) 
is compromised by this mismatch. 

On a more general point, full electromagnetic bound- 
ary conditions are not appropriate for MHD since the lack 
of displacement current means the MHD is pre-Maxwell; 
this should also be taken into account in any scattering 
description. 



III. CONCLUSIONS 

In this article we have illustrated some shortcomings 
of the magnetofluid modelling of the universe when used 
to determine magnetic field contributions to cosmic de- 
velopment. MHD necessarily can only describe slow 
timescale effects, simply because that is its mathemat- 
ical and physical basis. It is no surprise that MHD de- 
scriptions yield slow evolutionary behaviour; the model 
is incapable of delivering anything else. 

Assessing the influence of the magnetic field in post- 
decoupling cosmic structure based on such models is 
therefore at best approximate, and at worst very mis- 
leading. It is now timely to revisit the fundamental basis 
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of magnctogcncsis and evolution in order to extract the 
best possible interpretation of the observations. 

This article is designed to provoke discussion: we can- 
not fully exploit the true information content of cosmo- 
logical data if the community persists in attempting to 
model physical phenomena in an inappropriate mathe- 
matical framework. Some of these contentious assump- 



tions are highlighted here. We do not offer detailed al- 
ternatives in this article, since the construction of new 
modelling environments is a challenge for us all. 

LFAT acknowledges the financial support of the Lev- 
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for research funding. 
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